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Quantum error correcting codes have been developed to protect a quantum computer from de- 
coherence due to a noisy environment. In this paper, we present two methods for optimizing the 
physical implementation of such error correction schemes. First, we discuss an optimal quantum 
circuit implementation of the smallest error-correcting code (the three bit code). Quantum circuits 
are physically implemented by serial pulses, i.e. by switching on and off external parameters in the 
Hamiltonian one after another. In contrast to this, we introduce a new parallel switching method 
that allows faster gate operation by switching all external parameters simultaneously. These two 
methods are applied to electron spins in coupled quantum dots subject to a Heisenberg coupling 
H = J(t)S\ -S2 which can generate the universal quantum gate 'square-root-of-swap'. Using parallel 
pulses, the encoding for three-bit quantum error correction in a Heisenberg system can be acceler- 
ated by a factor of about two. We point out that parallel switching has potential applications for 
arbitrary quantum computer architectures. 



I. INTRODUCTION 

Quantum computers are capable of efficiently solving 
problems such as prime factoring |lj or simulating other 
quantum systems for which no efficient classical al- 
gorithm is known. A quantum computer is a device that 
stores and processes information which is physically rep- 
resented in its quantum state ||. Typically, such a de- 
vice contains a collection of quantum two-state systems, 
e.g. spin-1/2 particles. The state of each two-state sys- 
tem then represents a quantum bit, or qubit, the small- 
est indivisible unit of information in a quantum com- 
puter. Computations are driven by interactions between 
the qubits, generating logic gates operating on them. A 
quantum gate operating on M qubits can be represented 
as a 2 M x 2 M unitary matrix. Usually, a computation 
or algorithm is split up into a series of elementary gate 
operations involving only one or two qubits. In this repre- 
sentation, algorithms are also called quantum circuits. It 
has been demonstrated that there exist elementary two- 
qubit gates U which are universal when complemented 
with a sufficiently large set of single-qubit gates . This 
means that any quantum algorithm can be split up into a 
quantum circuit which contains only U and single-qubit 
gates. Quantum circuits are in general not the most ef- 
ficient way of implementing a quantum computation, as 
we will demonstrate in this paper. 

First experimental realizations of quantum computa- 
tion using trapped ions ||] , optical cavities || , and NMR 
0, involving two or three qubits, have been reported. In 
contrary to all of these systems, solid-state implementa- 
tions have the potential for a large-scale quantum com- 
puter involving hundreds or thousands of qubits. In this 
paper, we will concentrate on a theoretical proposal to 
use coupled semiconductor quantum dots in which the 



spin of the excess electron on each dot represents a qubit 
Q . Apart from electron spins in quantum dots, a number 
of other solid-state systems have been proposed for quan- 
tum computation: Nuclear spins of donor atoms in silicon 
[|L Josephson junctions p0| , d-wave Josephson junctions 
[ ll| , and charge degrees of freedom in quantum dots [fl2| . 
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FIG. 1. The circuit representation for three-bit quantum 
error correction, where time is evolving from the left to the 
right. First, the three qubits (represented by the horizontal 
lines) are initialized. The following unitary transformation U 
on qubit 1 prepares the state |$). The encoder E encodes the 
state |$) in an entangled state of all three qubits. In the next 
step, (simulated) decoherence partly disrupts the state. After 
the decoding step (which is identical to the encoding), the 
qubits 2 and 3 are measured. If they are both one, qubit 1 has 
to be flipped (C = a x ), otherwise qubit 1 is left unchanged 
(C — 1). If no more than one bit flip error occurred, the 
resulting state in register 1 is again |$) despite the presence 
of the decoherence. 

The physical implementation of quantum computation 
hinges upon the ability to find or design systems in which 
quantum phase coherence is maintained over long times 
compared to the duration of the typical controlled coher- 
ent operation. The discovery of quantum error correcting 
codes has been a landmark in the effort to find methods 
to protect a quantum computer from the decohering ef- 
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fects of a noisy environment |L3|] . The smallest quantum 
error correcting code for one qubit involving three code 
qubits has already been implemented in NMR E^] . 

In this paper, we present theoretical methods for find- 
ing an optimal implementation of three-bit error correc- 
tion. The optimization is understood here in terms of 
switching speed and switching complexity. The former is 
mainly motivated by the presence of decoherence which 
makes fast switching desirable, while the latter can be 
necessary if the physical implementation sets limits to 
the complexity of the switching. The two optimization 
goals usually are in conflict to each other, i.e. a fast im- 
plementation usually requires a complex switching mech- 
anism while switching with a simple mechanism is slow. 
First, we will study the "simple and slow" switching pro- 
vided by quantum circuits, and try to optimize it. Then, 
we will go on to "complex and fast" switching, for which 
we introduce parallel (as opposed to serial) pulses for the 
control parameters of the system, and show that the par- 
allel pulses allow faster switching than serial pulses. We 
also introduce a numerical method for finding such paral- 
lel pulses for arbitrary gates and Hamiltonians. We note 
that in a similar approach, Sanders et al. |l5| ] use a ge- 
netic algorithm to find optimized complex pulses to gen- 
erate quantum gates using an Ising type Hamiltonian for 
optically driven quantum dots coupled by dipole-dipole 
interactions. Our approach differs from this work in the 
underlying Hamiltonian (i.e. the mechanism proposed for 
quantum computation); in addition, the parallel pulses 
we suggest are general pulses that are discretized in time, 
whereas the pulses in jl5| are chirped Gaussian pulses. 

Three-qubit quantum error correction is able to cor- 
rect either one bit flip or one phase flip error (correct- 
ing both types of errors requires a code with at least 
five code qubits Q). The operation of the scheme is 
shown in the quantum circuit, Fig. [|. First, all three 
qubits (represented by horizontal lines) are initialized to 
the state |0), e.g. by polarizing the spins with a strong 
magnetic field. Then qubit one is rotated into an arbi- 
trary state |*) = U\Q) = a|0) The purpose of the 
quantum error correction scheme is to protect this state 
from external decoherence. In order to do this, the state 
|\&) is encoded in the highly entangled three-qubit state 

\9 L ) =^(|*)i|0) 3 |0) 3 )=a|0>i|0) a |0)3+/3|l)i|l)2|l)3. 

The encoded state can then be subject to exter- 
nal noise causing a (partial) bit flip of one of the spins, 
— > exp(ieaf — without the information 

contained in the encoded state being lost. The state 
is recovered (decoded) by applying the inverted en- 
coding network E^ 1 = E' = E which is identical to E 
since E is Hermitian. Then the qubits 2 and 3 are mea- 
sured (for physical implementations of quantum measure- 
ments on spins in quantum dots, see Refs. Pjlq-fLSfl. If 
both qubits are in state |1), then qubit 1 is flipped, oth- 
erwise it is left unchanged. This restores the state \^>) 
in qubit 1 which then can be measured in order to check 
the functionality of the scheme. 

For a first experiment one would probably want suffi- 



ciently low noise such that the state |^) is not destroyed 
by "natural" noise. One would then introduce bit flips 
"by hand" by applying a random oscillatory magnetic 
field in the x direction and then check whether those ar- 
tificial errors can be corrected for. 

For the sake of concreteness, we will apply our meth- 
ods to a system of coupled spins Si with s — 1/2 (each 
representing a qubit), subject to isotropic spin-spin in- 
teraction and local magnetic fields. With this model we 
capture the physics of electrons in coupled quantum dots 
H . We emphasize that a generalization to other systems 
with different Hamiltonians is straight-forward and does 
not require a new method for optimizing the switching 
process. 

Our paper is organized as follows: In Sec. [5], we in- 
troduce the formalism that we use to describe the dy- 
namics of electron spins in coupled quantum dots and 
other Heisenberg systems. The methods developed in 
Sees. Ill and IV, including the use of parallel pulses, are 
not special to the Heisenberg Hamiltonian Eq. (|l]). As 
an example we give some results for transversely coupled 
spins (XY model) in Section |V], because they are encoun- 
tered when electronic spins are coupled using cavity- QED 
p9| . The results of the Sections p^jv| are independent of 
the mechanisms that are involved in their physical imple- 
mentation - they are derived under the assumption that 
the model Hamiltonian Eq. (jl]) (or Eq. (ff6|)) is exact. 
In Section fv| we discuss some limitations and necessary 
cond ition s for the validity of this approach. Finally, in 
Sec. |VIlj , we give a detailed list of instructions for both 
serial and parallel switching which must be followed in 
order to implement three-qubit quantum error correction 
in a system of spins subject to Heisenberg interactions in 
experiment. 



II. MODEL 

In the system we consider the qubit is represented by 
the spin 1/2 state of the excess electron in a quantum 
dot, i.e. the "spin up" state | f) is identified with the 
logic state |0) = |t) and likewise |1) = where the 
quantization axis is chosen along the z axis, <J Z \\) = +|T) 
and a%) = - 1|>. 

The excess electron spins in a pair of quantum dots 
which are linked through a tunnel junction can be de- 
scribed by the Heisenberg Hamiltonian |qjL6| 



H(J, B.;, Bj) = J Sj • Sj + Bi • Si 



Bj • Sj , 



(1) 



where Si = Ci/2 describes the (excess) spin 1/2 on dot i 
and J denotes the exchange energy, i.e. the energy gap 
between the spin singlet and triplet states ||. This ef- 
fective Hamiltonian can be derived from a microscopic 
model for electrons in coupled quantum dots |Q, see 
also Sec. VI. It is found that J can be changed using a 
variety of external parameters. Tuning the gate voltage 
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between the coupled dots changes the height of the tun- 
neling barrier and therefore directly alters J. Note that J 
is exponentially sensitive to barrier changes. Also, apply- 
ing a magnetic field perpendicular to the 2DEG within 
which the quantum dots are defined greatly influences the 
exchange coupling J and can even result in a sign change 
of J for unscreened Coulomb interaction jl6| . Some cou- 
pling of the spin Si to a local external magnetic field B; 
is also necessary for quantum computation, and has been 
included in the Hamiltonian Eq. (jj). Note that we have 
included the factor gifJ-B in the definition of the mag- 
netic field Bi , where gi is the g- factor for dot i and /is is 
the Bohr magneton. The physical realization of the field 
gradients or inhomogeneous g-factors required for the lo- 
cal magnetic fields is challenging, but there exist several 
possibilities for generating them fll6| , [T8| |. 

From the Hamiltonian Eq. ([!]) we can generate the fol- 
lowing set of quantum gates, 



Ui(tf>) = exp(i</>- Sj) 
S = S(i,j 



U 1 ' 2 

w swap 



e-™' s exp 



Of 



s, ■ s. 



(2) 
(3) 



The single-qubit operation Ui{4>) for the spin Si is gen- 
erated by applying a magnetic field pulse B(t) at the 
location of the spin Si such that J drB(r) — <p. Simi- 
larly, the 'square-root of swap' gate [|| (which we denote 
by S in the following) is obtained by switching the in- 
teraction J(t) between the spins Si and Sj such that 

Jq cItJ(t) — 7r/2. We introduce the circuit notation for 
S in Fig. ||a. Note that Ui((p) is 47r-periodic in (ft, and 
27r-periodic up to a global phase —1, which for our pur- 
poses is not important. Eqs. (||) and (|J) are a univer- 
sal set of gates. Other powers U^ ap of the swap gate 

t/swap = S 2 = e-^exptwrSi ' s j)- (\ ab ) = W)i ® 1%) 
|00) i * |00), |01) ^ |10), |10) -> |01), |11) ^ |11), (4) 



can also be generated by the Hamiltonian Eq. (|l|) , but are 

7-1/2 
; swap 

1 /2 

Note that we cannot substitute J7 S wap by any other power 
of a 



not necessary for universality, once S 



swap 



because only £/^ ap with a half an odd integer 



has the ability to produce maximal entanglement from 
an unentangled state. 

We can use 'square-root of swap' to generate the con- 
trolled phase flip gate Ucpf- 



|00) i — ^ |00), |01) h-> |01), 1 10) h-> 1 10), |11) 



111), (5) 



with the quantum circuit depicted in Fig. pi (time evolv- 
ing from the left to the right), or formally |8|] , 

C/cpf = e^e^'e^^'Se^'S, (6) 
which in turn is related to XOR gate f7 X oR:(Fig. |b) 

|00) ^ |00), |01) ^ |01), |10) ^ |11), |11) ^ |10), (7) 
by the basis change 



Uxor = VU C p F V^ 
V = exp(-«7rSg/2) 



(8) 



(a) 



(b) 



FIG. 2. Circuit notation of two universal gates: (a) The 
'square-root-of-swap' (S) gate, (b) the XOR gate. 
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FIG. 3. A circuit representation for the conditional phase 
flip (CPF), Eq. (|), as given in Eq. (|). The single qubit 
rotations are F = e i7vS ' and G = e^ s \ The CPF is related 
to the XOR gate Eq. (Q) by the basis transformation Eq. (^) . 

Since XOR with one-bit gates is a universal quantum 
gate Q, this confirms that Eqs. (§) and (§) arc a univer- 
sal set of gates. In what follows, we will use the XOR 
gate to construct the gate E that performs the encod- 
ing for three-qubit quantum error correction, as shown 
in Fig. [l], and can be obtained by two successive XOR 
gates (Fig. |a), 



£ = C/xor(1,3)E/ X or(1,2). 
The very similar quantum gate (Fig. ^b) 
S t = C/ XO r(2,3)C/ X or(1,2), 



(9) 



(10) 



has the nice property that it can be used for implement- 
ing the quantum teleportation of one qubit as a quantum 
computation |^(|. It is clear that our analysis of the XOR 
gate can also be used for implementing this gate. 















J 



(a) 



(b) 



FIG. 4. The quantum circuits for (a) the three-bit en- 
coder E, cf. Eq. (Bf), and (b) the teleportation encoder Et, 
cf. Eq. ©. 
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III. SERIAL PULSE MODE 

In the foregoing discussion we made clear why it is 
desirable to generate certain quantum gates or networks 
such as XOR, E, and Et, and that it is indeed possible 
to produce them using a system of spins that are mutu- 
ally coupled by the Heisenberg interaction Eq. ([j]). In 
fact, we know that we can generate every quantum gate 
using those interactions, since Eqs. ^ and (||) explic- 
itly tell us how to produce XOR, which together with 
the set of single-qubit operations forms a universal set of 
gates for quantum logic . We now go one step further 
and ask ourselves which is the most efficient way of im- 
plementing a certain quantum gate. More precisely, we 
are interested in minimizing the switching time t s for the 
desired quantum gate. This kind of optimization is cru- 
cial because the error probability per gate operation is 
proportional to the switching time, e = t^/t^, where 
denotes the dephasing time of the system. Other criteria 
for optimization can be added, if e.g. one kind of elemen- 
tary task (say, spin-spin interactions) is much harder to 
perform than another (such as single-qubit rotations), or 
if the switching of parameters turns out to be difficult. 

For the purpose of finding an optimal implementa- 
tion of quantum gates, we first define which set of el- 
ementary operations we are going to use. We will call 
this set the serial pulse operations, since they can be 
achieved by 'switching on' exactly one of the parameters 
p = (J, B!,B 2 ) in the Hamiltonian Eq. (|l|) for some fi- 
nite time. Clearly, XOR does not belong to this class 
of gate operations - it takes the whole sequence Eq. (||) 
to produce it. Note that the definition of serial (and 
later parallel) pulse operations depends on the Hamilto- 
nian and on how we parametrize it. The distinction only 
makes sense if serial pulse operations correspond to phys- 
ically switching on and off a part of the device, e.g. the 
magnetic field at the location of one of the spins. We will 
use the serial pulse operations defined in Eqs. (0) and (0). 



A. XOR gate 

As a first example for our efficiency analysis we take the 
sequence Eq. (||) for XOR. We do not try to optimize the 
length of the single-qubit pulses. Instead we investigate 
whether it really takes two instances of S or if XOR can 
be performed with one S plus single-qubit operations. 
This is most reasonable if coupling two qubits is more 
costly (e.g. due to decoherence) than operating on a 
single qubit. Thus, our question is whether 



Si = [u2i <8> u 2 2\ S [un <g> W12] 



(11) 



is equal to Uxor f° r some choice of single-qubit gates 
u nm = U m ((f> nm ) or not. We will shortly prove that the 
answer is negative and it indeed takes at least two S to 
produce XOR, but first we introduce the method we de- 
veloped in order to prove this kind of 'no-go' theorem. 



The idea is that very often quantum gates can be distin- 
guished by their ability to produce entanglement. This 
property of quantum gates has the advantage that it is 
invariant under concatenation with arbitrary single-qubit 
gates. 

We denote the product (pure) states in our 
two-qubit Hilbert space U = Hf 2 by V = 



{l*> 



G H 



W - \<P) ® \x)l \<p), \x) € U 2 ) [[§. Here, H 2 
denotes the single-qubit Hilbert space with basis |0), |1). 
A state \%>) ^ V is called entangled. For every quantum 
gate (unitary matrix) U acting on TL 1 we define the sub- 
set V(U) = e V U\i&) e 7>j C V of product states 

which are mapped onto a product state by U. The idea 
is now simply that two quantum gates U\ and U% which 
have different P-sets, V{U\) ^ TiJJ-i), obviously must be 
different: U\ 7^ U2 (note that this implication cannot be 
reversed). The V-set of XOR is 



V(U XO r) = {\0cf>),\l<p},\<i>±)\\<t>) ew 2 } 



(12) 



where we used the notation |±) = (|0) ± |1})/V2. 

In order to find V(S) it is useful to convince oneself 
that 'square-root of swap' operates on a product state 
\4>X) according to the very intuitive formula 



1 



S\<f>x) = T - -.{\4>x)+i\x<t>)) 
1 + 1 



(13) 



by first checking it for the basis of products of |0) and 
|1), and then using that the right-hand side of Eq. ( |l3| ) 
is linear in \<p) and |x). From this rule we conclude that 
all product states become entangled by S unless they are 
the product of two equal single-qubit states, 



T(S) = {\ 



= 1 1*) eV 



G H 2 



m = 



(14) 



For any choice of the u nm in Eq. (|11|) , we can construct 
the state |0) ® u| 2 wn|l) which is in V(Uxok) but not 
in V(Si) since is entangled. Therefore, V{Si) ^ 

V(Uxoa) and consequently Si ^ Uxor for any choice of 
Unm- Thus, the sequence given in Eq. (|^) is optimal in 
the sense that both 'square-root of swap' operations are 
really needed. Allowing arbitrary powers of J7 swa p does 
not reduce the number of two-qubit gates required for 
XOR either, since one USL, a where a is not an even 

1 swap " 

multiple of 1/2, cannot act as a perfect entangler which 
is required for the XOR gate. For completeness, we give 
here the generalization of Eqs. (|l^) and ( |l4] ) for arbitrary 
powers of swap, 



e -i7ra/2 CQS 



— j \4>x) +iem 



V(S), a ^ integer, 
V, a = integer. 



1X0)) , 
(15) 
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B. Three bit encoder E 



U 2 = [u® V(2,3)]S(1,3)E/'W, 



(19) 



Regarding the three bit encoder E, Eq. (||), our result 
tells us that the straight-forward implementation of E 
requires 'square-root of swap' four times, i.e. twice for 
every XOR. This does not mean that there cannot be a 
more efficient implementation of E than given in Eq. 
We can try to implement E using the serial pulse gate 
set instead of XOR's. It turns out that this is impossible 
with fewer than four S gates. The analysis still relies 
on the previously introduced V-set but is slightly more 
complicated than the one for XOR since in the case of 
three qubits each gate S can be applied to one of three 
possible pairs of qubits. 

It is clear that just one use of S (plus arbitrary single- 
qubit operations) cannot produce E, 

Ui = [u2i <8) u 2 2 <8> U2a]S(i, j)[un ® u i2 <g> W13] ^ E, (16) 

for any choice of u nrn for the simple reason that E is 
able to entangle the qubit 1 with 2, and also 1 with 3, 
whereas S(i,j) can only entangle the qubits i and j with 
each other (at most one pair). 

It is less obvious that none of the sequences 



U 2 = U {3) S(k,l)U^S(i,j)U {1) 



(17) 



with = u n i (g) it„ 2 &> Un3 can reproduce E. The 

idea of the following argument is the same as for the 
one for XOR: We are seeking a state |\E') that becomes 
entangled when acted on with the operator U2 given in 
Eq. ( |17| ) but remains unentangled under the operation 
E, i.e. I*) G V(E) \ V{U 2 ). This is the case if |*> is 
entangled by S(i,j) and not disentangled by S(k,l). We 
can exclude the case where (fc, I) = (i, j) or (fc, I) — (j, i), 
using the same argument as for U\, defined in Eq. (16). 
In the remaining cases it is clear that if |^) G V(E) is 
chosen such that it is entangled by S(i,j), then it will 
not be disentangled again by S(k,l), and we are done. 
Since Eq. ( p"7j ) is invariant when i and j (or k and I) 
are interchanged, we can always arrange that i, j, and 
k are mutually different, and I — i. In the case where 
j ^ 1 the state |*> = |0), <g> u^tn^l^ <g> |0) fe G V{E) is 
not in T > (U 2 ), because the entanglement between qubit 
j and the qubits i and k created by S(i,j) cannot be 
undone by S(i, k). For j = 1 we choose the state = 
M^Mijll)^ ® |0)j ® |0)fe G "P(E) with the same property. 
This concludes our proof that there is no circuit U2 of 
the form Eq. (|l7]) which is equal to E. Note that this 
conclusion is independent of the choice of single-qubit 
operations in U2, hence the inequality we proved concerns 
all circuits of the type U 2 - 

As an example consider the circuit 



U 2 - t/ (3) S(2,3)[/ (2) S*(l,3)[/ (1) , 



(18) 



for which i = 3, j = 1, and k 
the form shown in Fig. ^, 



2. By rewriting U2 in 



with?; = -U31M21 and F(2,3) = (1<E>U32®U3 3 )S(2, 3)(1® 
M22 ( 8'W23), and using Eq. (13h, we observe that the prod- 
uct state I*) = |0)i®|0) 2 ®Mi 3 Wn|l)3 G V(E) is mapped 
to 



1 



1 + * 



w|a>i ® V|tf) 23 + ® ^17)23 ■ (20) 



The unitarity of un implies that |a) = wn|0) is orthog- 
onal to = wn|l), and I7) = uia|0) ® |a) is orthogonal 
to I <5) = U12IO) (8> \ f3). The gates u and V are also unitary, 
thus v\a) is orthogonal to and VW) is orthogonal to 
V\S), which implies that L^l^), Eq. (p0|), is an entangled 
state between subsystem 1 and subsystem 2 and 3. From 
this we conclude that |\&) G V(E) is not in T'(U2), and 
therefore E U 2 - 



u n 



U12 



U13 



U21 U31 



V 



U22 -f — U32 



V 



FIG. 5. Quantum circuits of the type described in 
Eq. Dashed lines represent the grouping in Eq. (|l9[). 

Next, we develop a proof that even with the use of 
three S gates, E cannot be implemented. Since each S 
can couple one of three possible pairs ik,jk — 1 , . . . , 3, 
ik 7^ jk, of qubits, there are 3 3 = 27 sequences including 
three 'square-root of swap' (S) gates, having the form 

U 3 = uWS(hj3)UWSfa,h)& 2) S(hJi)uW, (21) 

with arbitrary single-qubit gates = u n \ ® u n 2 ® 

u n z- First we observe that if (12,^2) = (*3jj3) or 
(*2>J2) = (isj^a) then we can apply the same argument 
as for circuits of the type U2 with = (*i, ji) and 

(k,l) = (12,33)- In the case where the first two S gates 
(but not the third one) act on the same pair of qubits, 
(kji) = (h,32) or = (j2,fa), we note that ei- 

ther |0) il |0) J - 1 |0) fe G V(E) or |0) il |0) il |l) fc G V{E), where 
k ii,ji, becomes entangled by U3. Therefore U3 7^ E if 
the first two or the last two S gates operate on the same 
pair of qubits. In all other cases, we can label the three 
qubits with three distinct numbers a, b, and c between 
1 and 3 such that S(ii,j±) = S(a,b), 5 r («2,j2) = S(a,c), 
and S(i3, j'3) = S(b,x), with x = a or x = c. The state 
I*), defined as 
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|0) o <8> u\ b Ula\0)b <8> u\ c U2 c U2aUla\l)c, (if a = 1), 

u\a U lb\0)a <8> |0)fc ® uL4c U 2aWlh|l)c, (if & = 1), 

M la M 2a U 2cWlc|0)a ® uJ b 4a U 2 C -"l c |0) 6 <g> |l) c , (if C = 1), 

(22) 

is chosen such that G 'P(-E') and has the 

property that 5(a, 6)f/W|$) is unentangled, but in 
S^a, c)U^S(a, b)U^\^) there is entanglement between 
the qubits a and c. Finally, the S(b, x) gate cannot re- 
move the entanglement; in the final state U-^), either 
the qubit a (if x — c) or the qubit c (if x = a) is entan- 
gled with the other two qubits, thus ^ V(JJz). Since 
|*) G V(E), this concludes our proof of the statement 
[^3 7^ _E. In order to illustrate our proof, we apply it to 
the specific example (a = 2, b = 1, c = 3, a; = c = 3) 

J7 3 = £/ ( %(l, 3)?7 (3) 5(2, 3)[/ (2) S(l, 2)f/ (1) , (23) 

which can be written in the form (shown in Fig. ^ 

WS(2,3)U (2) S(1,2)U W , (24) 

where W is a gate that does not couple qubit 2 with any 
other qubit. Applying this operator to the state = 
1 0> i <g> u\ 2 un\0)2 ® u\ 3 ul 3 U22Un\l)3 proves that U 3 ^ E 
for any choice of the U^ n \ because U^\^) is entangled, 

I*) i v(u 3 ). 



cannot entangle any pair of the three qubits, as it should 
in order to reproduce Et- 

A circuit of the form U2, Eq. (|l7|) , involving two S 
gates cannot be equal to Et either. We first note that if 
S(i,j) — S(k,l), then U 2 cannot produce entanglement 
between the third qubit m 7^ i,j and i or j, and thus 
U2 =/= Et, because Et can entangle any pair of qubits. If 
S(i,j) £ S(k,l), we find that either |*) = |0) l |0) J |0) fe G 
V{E T ) or I*) = |l) l |0) J |0) fc G V(E T ), where k + 
is entangled by the gate S(i,j). Since this entanglement 
cannot be undone by S(k,l) 7^ S(i,j), the state L^l^) 
is entangled. Thus, we have found a state which is in 
T(E T ) and not in T(U 2 ) and therefore U 2 ^ E T . 

We finally explore whether a circuit U3 containing 
three 5 gates as in Eq. (HJ) can reproduce Et- For 
S (12,32) = S(i,3,j3) we can see that this is not the case 
by applying the same argument as above for a circuit 
with two S. In the opposite case, S(i2,32) 7^ S(is,jz), 
either the state 10)^0)^ |0) fc G V(E T ) or the state 
10)^10)^11)*: G V(Et), with k f is entangled by 

S (12,32) or S(i 3 , j 3 ). Because S(i 3 , j 3 ) ^ S(i 2 J 2 ), entan- 
glement produced by S(i 2 ,j 2 ) is not undone by S(i 3 ,j 3 ) 
and therefore we have found a state in V(Et) which is 
not in V{Uz), completing the proof for U3 7^ Et for the 
case of three S gates. This finally implies that, like E, 
the teleportation encoder Et cannot be constructed us- 
ing fewer than four S gates. 



2 U 12 



U21 



3 U 13 



^23 



U32 



U 41 



U42 



W 



FIG. 6. Quantum circuits of the type described in 
Eq. ([23|). Dashed lines represent the grouping in Eq. ([m]). 

C. Teleportation encoder Et 

Now we consider the teleportation encoder gate Et, 
Eq. jl0|), which is shown in Fig. [IJd. The gate Et consists 
of two XOR gates, but in contrast to E these two XORs 
are arranged in a less symmetric way. As a consequence, 
V(Et) is only a subset of (and not equal to) V(E), since 
e.g. |0 + 0) is in V(E) but not in V(E T )- Again, we can 
ask whether it is feasible to assemble Et using less than 
the four S gates that are used when we simply combine 
Eqs. (||), (g), and ([To|). The answer is again negative. 
The proof of this statement is similar to the one given 
for E. It is again clear that a circuit involving one S 



D. Numerical search 

The method that was presented for proving inequal- 
ities between two gates U\ and U2 involving the sets 
V(U\) and V(U2) has the advantage that it yields rigor- 
ous results although we do not know the details about 
the involved single-qubit operations. Sometimes how- 
ever, proofs become lengthy and rather unsystematic, so 
we would like to have a better tool for the complex cases. 
Unfortunately, we do not have such a tool which is ca- 
pable of giving rigorous proofs for inequalities, like the 
P-set method. However, we have developed a computer 
algorithm that searches for the M qubit gate U g in the 
set 

[/{[/(»)} = uW +r >X N uW ■ ■ ■ X 2 U^X 1 U { - 1 \ (25) 

where the X n are arbitrary but fixed 2,3, ...,M qubit 
gates and U^ 71 ' — u n \ ® ■ ■ ■ <g> u n M ar e arbitrary and 
variable products of single-qubit gates (several u n k can 
be unity). A result of a numerical search is a list 
{link} of single-qubit gates, which satisfy the equation 
U{W n >} = U g . The computer algorithm can therefore 
(in the case of a successful run) 'prove' equalities, but 
one or several unsuccessful runs do not constitute a proof 
that the gate U g cannot be constructed with a given se- 
quence X n) n = 1, N. Note that the situation is thus 
exactly opposite to the P-set method. The operation of 



G 



the computer algorithm consists of minimizing the func- 
tion 



f({u nk },a) = \\e la U{u nk }-U g \\ 



(26) 



numerically in the space of all possible combinations of 
single-qubit gates u n /-, where the matrix norm is given 
by ||A|| 2 = Tr(A^A). The single-qubit gates u nk are 
parametrized by the three angles tp n k, 9 n k, and ip n k ac- 
cording to the prescription 



cos(6» nfc ) 
e^»* sin(0„ fe ) 



aiiVnk+lpnk) 



iip„ 



cos{V nk ) 



(27) 



Including the global phase a we count 3M(N+ 1) + 1 real 
parameters. If the numerical search yields a minimum 



f({u nk },a) = 0, 



(28) 



then the corresponding sequence Eq. (129) is identical to 



The numerical results for U g = Uxor and X n 



S 



support the analytical result, i.e. there is no circuit for 
N = 1. For N — 2, we find a vast number of circuits other 
than Eq. (||) combined with Eq. (^). For U g = E, Et, we 
do not find a solution for N < 4, as guaranteed by the 
result of our previous analysis. 

The impossibility of reducing the number of S gates 
required for E led us to the idea that it might be 
useful to replace S by a three-qubit gate which is di- 
rectly generated by the three-qubit Hamiltonian #3 = 
J(Si • S2 + S2 • S3 + S3 • Si), describing three simultane- 
ously interacting spins with equal coupling constant J. 
We find that the analogue of S for three spins is the gate 



5 3 = e -/ 3 exp(^), 



(29) 



which is obtained when the interaction J is switched on 
for time r = 4?i7r/9J. In analogy to Eq. ( |l3| ) we can 
express the action of S3 on product states as 



creS 3 



(30) 



where the second term is the symmetrization of the in- 
put state (1S3 denotes the permutation group of three ob- 
jects). Whereas S 4 = 1, we find that S3 = — 1. In exact 
analogy to the 'square-root of swap' gate, V(Ss) consists 
of states that have the form |a) ® |a) ® |a), \a) S H 2 - 
We can show that E is not equal to any gate involving 
only one or two S3 . The case of three S3 was studied nu- 
merically, but no circuit representation for E was found. 



IV. PARALLEL PULSE MODE 

Operating a system described by the Hamiltonian 
H(p) with parameters p = (J, Bi,B 2 ) given in Eq. (Ill) 



as a quantum gate in the serial pulse mode is not opti- 
mal in the following sense: If several or all parameters p 
can be changed simultaneously, we expect that a given 
quantum gate, say XOR, can be performed faster than 
by changing only one parameter at a time as in the serial 
pulse mode. Generally, all parameters p are arbitrary 
functions of time such that the time evolution operator 
after time t is a functional in p given by the time-ordered 
exponential 



u t W)\ 



TexpQ jf H(p(T))dr 



(31) 



Given some quantum gate U g , we would now like to solve 
the integral equation Ut[p{r)] — U g for the functions p(t) . 
For unrestricted time t and unbounded functions p(r) we 
immediately know how to construct such a solution by 
using the known universal set of gates Eqs. (||) and (||) in 
the serial pulse mode. In general, this is not the optimal 
solution of Ut\p[r)] — U g . An optimal solution is given 
by a set of bounded functions |pi(r)| < Mj requiring 
minimal time t for a fixed set of bounds Mi. Since it is 
not feasible to find an optimal solution among all such 
bounded functions, we will restrict ourselves to piecewise- 
constant functions. Splitting up the time interval t into 
N > 1 parts, we write 



U N {p 



(i) 



^U N (p^)---U 2 (p^)U 1 (p^), 



Uk(p {k) ) =exp{2^i?(p (fc) )} 



(32) 



S 3 |aia 2 Q!3) = -e tn/3 \a!a 2 a 3 ) + V] \a al a a2 a<j?,) , the relation 

2V3 



For every time "slice" , we have the freedom to choose a 
new set of parameters p^ — (J, B 1 ,B 2 ). Note that we 
allow for an arbitrary total phase </>. By discretizing the 
problem in this way we have reduced the free parame- 
ters in the problem from the P functions Pi(t) to PN 

(k) 

real parameters p\ , i = 1, ...,P, k = l,...,N where P 
denotes the number of parameters p (in the case of the 
Heisenberg Hamiltonian Eq. (Eh, P = 7). The functions 

Pi(t) are related to the discrete parameters through 



Pi{t) 



2nh 

Tk 



-P 



(fc) 



tk-1 <t<tk, 



(33) 



where Tk — tk — ife-i and to — 0, — t; the time step 
Tk has been absorbed into the dimensionless parameters 

(k) 

p\ . Once the problem is discretized, it becomes suitable 
for numerical treatment using the minimizer algorithm 
presented in Section III, minimizing the function \\U g — 

Ulf{Pi , 0}|| 2 with respect to the PN+1 parameters p\ k ^ 
and (f>. We try to find a solution to J7jv{p$ ,<^>} = U g 
starting from N = 1 and then increasing N in unit steps. 
In practice, N is limited by the available computational 
resources. 

One approach to the problem would then be to fix 
N and U (e.g. use time steps of equal size, Tk = 
t = t/N). Then, the constraint |pi (x) | < Mi implies 
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\Pi°^\ < T kMi/2Trh. In the following, however, we solve 
UnIp^ 4>} — U g with fixed N (chosen as small as pos- 
sible) without any constraint for p^' and then calculate 



„7riA 



-TtiA 



t for given bounds M, using the formula (cf. Eq. (|33|)) 



N 
fc=l 



N 



2irh 

m7 



(34) 



The parameter p^ k ' with the largest p\ k ^ /M^ ratio deter- 
mines the switching time for the fc th step. 



A. XOR gate 

We now want to use this method for finding a pulse 
sequence p 1 -^ that generates the quantum XOR gate, 
Eq. (Q). Since XOR is the same as the conditional phase 
flip (CPF) up to the basis change Eq. (g), we will first 
try to generate CPF. In the S z basis, the Heisenberg 
Hamiltonian H(J, B x , B 2 ), Eq. (Q), can be written in the 
following matrix form: 



Bf + Bf Bf - %B\ Bf - iB\ 

II Bf+iBf Bf-BI-J J Bf-iBf 

2 ( Bf + iB\ J Bl -Bf -J B. 

Bf + %B\ B. 



iB\ 



%B\ -Bl - Bl 



(35) 

where an irrelevant constant energy contribution is omit- 
ted. We find analytically that CPF can be obtained in 
one time-step (N — I), i.e. for constant parameters p, 

C/cpf = exp [2mH(J,B x ,B 2 )], 

J = k — n — 2m , 

2 

Bi = ^(0, 0, n + i + -\/fc 2 — J 2 ), (36) 
B 2 = i(0,0,n+i- Vk 2 - J 2 ), 

where n, and m are arbitrary integers, and k is an integer 
satisfying 2\k\ > \n + 2m + h\. This solution is obtained 
by setting Bf = Bf = 0, and diagonalizing the resulting 
Hamiltonian matrix 



H 





( A 















-J + B 


J 






1 


J 


-J-B 


1) 




V o 











(37) 



where A — Bf + Bl and B = Bf — B%. The conditional 
phase flip CPF is invariant under the basis change that 
diagonalizes Eq. (|37|), and we can then solve the equation 
exp(i^) exp(2niH) — C/cpf in the basis where both H 
and C/cpf are diagonal. This yields the four equations 



e where A± 



(-J± VJ 2 + B 2 )/2 and ±A/2 are the eigenvalues of H. 
From these four equations we obtain the result Eq. ( |3^ ) 
for 0, A = Bf + Bl, B = Bf - Bf , and J. Applying the 
basis change V from Eq. (ph to these solutions, we can 
build XOR with a total of N = 3 steps. The integers 
fcj m, and n can be chosen such that the switching time 
Eq. (HH) is minimal for a given set of constraints Mj, 
Mb x , Mb 2 - In the specific case where all constraints 
are equal to M, we find that the solution for k = 1, 
m = n = 0, 

J = \, Bf = ±(l + V3), B| = i(l-V3), (38) 



has the shortest switching time, 

^ = ^(1 + ^=0.683^, 



M 



(39) 



less than half the time which is used for the serial pulse 
quantum circuit Eq. (||), icPF,s = 1-5 • 2irh/M. Note 
that since the coupling is isotropic, the same gate (in a 
rotated basis) can be achieved with a magnetic field along 
any desired direction. In order to obtain the XOR gate, 
we must spend in addition twice the time 0.25 • 2Trh/M 
for the basis change V, Eq. (|§|). The total switching time 
is then 

2nh , ^ 2nh 
W P = ^(3 + \/3) = 1-183—, (40) 

about 59% of the time required for the serial pulse quan- 
tum circuit, Eq. (^), including the change of basis Eq. (||), 
ixOR.s = 2 • 2-KhfM. Of course, the basis change applied 
here is again a 'serial pulse' action and therefore not op- 
timal. We therefore study Eq. ([32]) directly for the XOR 
gate, without using the CPF gate. It turns out that no 
solution exists for N = 1. For N = 2 our optimizer 
algorithm finds the numerical solution 



Uxor 



exp 2niH(p 



(2)> 



exp 2TtiH(p 



(41) 



with the parameter values 



k 


j{k) 




B w 










1 


0.187 


-0.025 


0.464 


0.205 


0.195 


-0.420 


0.395 


2 


0.617 


-0.220 


0.345 


-0.384 


0.244 


0.353 


0.108 



(42) 

and the global phase 4> = —0.8481 • it. The total switch- 
ing time for equal bounds is in this case ixOR,p = 
(0.4643 + 0.6170)27T?i/^ = 1-0813 • 2nh/M, compared to 
tx(m,s = 2-2-kTl/M for the serial switching. The numbers 
in boldface in Eq. (^) indicate which parameter limits 
the switching time in each step. The solution Eq. (|4T 
appears to be a unique optimum for the case N = 2. 
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B. Three bit encoder E 

We can further parallelize the three-bit encoder E, 
Eq. (^). Instead of concatenating two XOR gates (which 
may or may not be produced using parallel pulses) we 
now try to find a more efficient parallel pulse sequence 
for E, given a system of three qubits which exhibit pair- 
wise couplings among each other that can all be switched 
on simultaneously. The Hamiltonian for this three-spin 
system can be written as 



H = 



E 

l<i<j<3 



^ B, • S,; 



(43) 



We find that there is a representation of the three-bit 
error correction encoder E which consists of three parallel 
pulses only, instead of the four which it takes to perform 
two sequential XOR gates. The following parallel pulse 
sequence produces E up to a global phase cj) = it/ 2: 



k 


T {k) 
J 12 


23 


J 13 


B lx 


d(«) 


B 3x 


1 


0.0000 


8.2500 


0.0000 


1.1153 


6.1737 


6.1739 


2 


-0.9256 


-5.3608 


0.7863 


5.7603 


5.2422 


1.7475 


3 


0.0000 


-1.7500 


0.0000 


0.4345 


3.5255 


3.5255 


k 










B W 

2z 


3z 


1 


-1.6737 


-0.2263 


0.2262 


1.1153 


1.4649 


-1.4649 


2 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


3 


2.1709 


1.4118 


1.4118 


0.4345 


1.2560 


1.2560 



(44) 

For equal bounds the total switching time of tg p = 17.54- 
2nh/M is much larger than the 4-pulse time 2ixoR,p = 
2.163 • 2-kTl/M . Note that a better 3-pulse solution was 
not found, but cannot be excluded. 



V. ANISOTROPIC SYSTEMS 

Systems where the spin-spin coupling is anisotropic are 
not described by the Heisenberg Hamiltonian Eq. (fy) that 
we studied as a generator for quantum gates in the previ- 
ous sections. In the two most notable cases, the Ising and 
the XY systems, it is known that universal quantum com- 
putation is possible. In the case of a system described by 
the Ising Hamiltonian H] = JSfS^ and a homogeneous 
magnetic field in z direction, there is a particularly sim- 
ple realization of the CPF gate with constant parame- 
ters, namely £7 C pf = exp(i7r(l - 2S{ - 2Sf + iSf 5|)/4) 
||. One might be tempted by this to 'transform' the 
Heisenberg interaction Eq. ([!]) into an Ising interaction 
by adding time-dependent fields Ho(t) — Bi(£) • Si + 
B 2 (i) • S2 to the coupling Hamiltonian V(t) — J(i)Si • S 2 
such that the coupling in the interaction picture, Vi (t) = 
U{t)V{t)U{t)\ with U{t) = Texp(i /*fr (r)dr), would 
be identical to the Ising coupling Vi(t) = Hi, or even 
to switch the coupling off and on using this method, i.e. 



Vi(t) — for a certain choice of Ho(t). It turns out that 
this is impossible, since the 'transformed' coupling must 
have the form 



Vi{t) = J(t)Si ■ (R(t)S 2 ), 



(45) 



where R(t) is a time-dependent rotation matrix. For J ^ 
0, this clearly excludes the complete 'switchin g off ' of the 
interaction. Furthermore, the coupling Eq. ( [45| ) is still 
isotropic at every instant t. 

In spite of the impossibility of using a Heisenberg sys- 
tem as an effective XY (or Ising) system by adding time- 
dependent fields, there are XY systems in nature which 
have been proposed for quantum computation pi| . In 
the case of Ising systems we have seen above that there is 
a very simple prescription for generating the XOR gate. 
We devote the rest of this Section to demonstrating that 
XOR can also be obtained with XY coupling. For two 
spins with s = 1/2 the XY Hamiltonian is given by 



Hxy = J(S'tS x 2 + S'fS v 2 ) = ~ 



/oooo 

10 
10 

Vo 



(46) 



where for the matrix representation we chose the S z ba- 
sis. The corresponding time evolution operator is 



1 



E/xy(</> = Jt) = exp(itH X Y) 



(47) 



1 



There is a qualitative difference between two qubits cou- 
pled via an XY and ones coupled by a Heisenberg in- 
teraction: it is impossible to generate powers f/^ ap 
(0 < a < 1) of the swap gate Eq. (||) with only one 
use of Uxy{4>) together with single-qubit operations. In 

particular, this is impossible for the 'square-root of swap' 
1/2 

gate E/swap- In spite of this, we found that the CPF gate 
can be produced by the serial-pulse sequence 



_ »tt/4 2i7rm-Si/3 2i7m 2 -S 2 /3 

x C/ xy (7r/2)e" s i[/xy(vr/2)e" 



-virSf/2-iirS$/2 

C- 1 



(48) 



where n x = (1, -1, l)/VZ and n 2 = (1, 1, -1)/V3. This 
is a proof that XY systems with single-qubit interactions 
are in principle capable of universal quantum computa- 
tion. 

We now consider parallel switching with the XY dy- 
namics, 



Hxy.b = Hxy + Bi • Si 

UxY,B(t) = exp(itH X Y,B)- 



B 2 • S2 



(49) 
(50) 



As in the case of Heisenberg interactions, we first con- 
sider the CPF gate which can be used to assemble the 
XOR gate as shown in Eq. (||). We have not found a pos- 



sibility to generate the CPF gate Eq. ( 



Hamiltonian with applied magnetic fields with constant 



J) with the XY 
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parameters (N — 1) using a numerical search [Q. If the 
switching is performed in two steps (N — 2), we find nu- 
merically that there are several possibilities to generate 
Ucpf hi the form 



U C p F = e l '»U 2 U 1 , 
2 m H XY . B {J {k \B^\B^; 



(51) 



1,2. 



where Uk = cxp 
Note that all magnetic fields can be chosen homogeneous 
(B^ = = B' fc )) and perpendicular to the y-axis 

(By = 0). Here we give one possible realization which is 
found numerically (4> = — 3ir/4): 



k 


j{k) 




B (k) 


1 


0.7500 


0.7906 


0.5728 


2 


0.5000 


0.0000 


0.2500 



(52) 



The total switching time for CPF, assuming equal 
bounds Mj = M B = M for J and B, is t^ Fp = 
1.291 • 2tt%/M, compared to t^ F s = 2.167 • 2irh/M for 
the serial pulse sequence defined in Eq. (jig) . 

In order to produce the XOR gate Eq. (0) we can im- 
plement the basis change Eq. (||) using the single-qubit 
rotation V. This procedure requires a total of four steps 
for the XOR gate. Another way of achieving XOR is 
the following sequence which we found numerically and 
which takes only three steps: 



U X OR = exp@iv/4)U 3 U2U 1 , 
with the following parameters: 



(53) 



k 


j(fe) 


B\ x 


B 2x 


B iy 


D* 

B 2y 


B\z 


BjS 


1 


1.802 


0.615 


2.045 


0.020 


0.316 


0.794 


0.130 


2 


3.344 


0.348 


0.718 


0.259 


0.493 


1.583 


1.062 


3 


1.903 


1.193 


0.705 


0.413 


-0.305 


0.589 


0.604 



(54) 

The total switching time of t x g R = 7.29 • 2nh/M (com- 
pared to 2.67 • 2-Kh/M using CPF and a basis change) 
indicates that Eq. (|54|) is not an optimal solution. 



VI. REQUIREMENTS FOR PARALLEL 
SWITCHING 



The parallel switching mechanism presented in the Sec- 
tions |^ and ^ relies on the following essential assump- 
tions: 

(a) Each of the parameters in the Hamiltonian can be 
varied independently. That is, the coupling can be varied 
independent of the magnetic fields in the Hamiltonian 
Eq. (0). 

(b) We know the exact relation between the externally 
controlled parameters (such as the electric and magnetic 
field or a gate voltage) and the parameters in the Hamil- 
tonian. 



(c) The switching is synchronous, with all parameters 
Pi varying with the same time profile Pi(t) = pif(t). The 
change of parameters does not have to be step-like, but 
can be chosen to have some smooth pulse form. Also, 
any pulse magnitudes pi are allowed. 

Whether the above requirements can be fulfilled de- 
pends on the underlying microscopic mechanisms which 
are responsible for the effective Hamiltonian, such as the 
Heisenberg Hamiltonian, Eq. (|l]). In our previous work 
[flq we have used the model Hamiltonian 



H 



2m ^-^ 

»=1,2 



-A(r 4 )) 2 +exiE 
c ) 



4a 2 



Vi 



« ri - r 2 



(55) 



with A(r) = B(—y, x, 0) /2 to describe the orbital dynam- 
ics of electrons in coupled quantum dots. Here, and 
Pi denote the location and momentum of the electron i 
which is moving in two dimensions in a double-well po- 
tential V with characteristic energy Tiloq and a magnetic 
field B perpendicular to the 2D electron system and an 
electric field E parallel to the coupling axis of the two 
wells. The distance between the quantum dots is denoted 
by 2a, the effective mass and the charge of the electron by 
to and e, and the dielectric constant of the material by k. 
As we pointed out earlier p6[ , the spin-orbit interaction 
H so = (luq i '2mQC 2 )S ■ L is very small for an electron in a 
parabolically confined quantum dot. Note however that 
this expression for the spin-orbit coupling contains the 
bare electron mass mo instead of the effective electron 
mass to, making the spin-orbit coupling in GaAs smaller 
than estimated in |l6| by a factor of Too/to ~ 15 p5[ . 
For a quantum dot with confining energy Tiloq = 3 meV, 
we obtain H so /Tiujq w 10~ 8 . 

Concerning condition (a) we have found that the spin- 
spin coupling J can be controlled by several external 
"knobs". The gate voltage V applied between the cou- 
pled quantum dots controls the height of the barrier for 
tunneling of an electron from one dot into the other and 
therefore strongly influences the exchange coupling J be- 
tween the electronic spins. In a similar manner, J de- 
pends on the inter-dot distance 2a. We have also found 
p6| that an external magnetic field B perpendicular to 
the 2DEG causes a strong change (even a sign reversal) 
of J. Not surprisingly, an electric field E applied along 
the coupling direction of the dots also changes the ex- 
change coupling, which can be understood as an effect of 
level detuning. When switching on a magnetic field, the 
effect of the field on J could be compensated by changing 
another independent control parameter, e.g. the electric 
field. In practice, one has to know the functional depen- 
dence J(V, a, B, E) in the range where it is used, see also 
(b). 

While a magnetic field perpendicular to the 2DEG 
strongly influences the exchange J, we can argue that 
sufficiently weak in-plane magnetic fields have little in- 
fluence on J. Classically, the motion of a particle in a 
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plane is not affected by a magnetic field in the plane, since 
the Lorentz force is orthogonal to the plane. Quantum 
mechanically, we can describe a particle in a magnetic 
field confined to a plane by the Hamiltonian 



H 



1 

2m 



(56) 



where the vector potential A = B(Q, — z, 0) corresponds 
to a magnetic field of magnitude B along the x axis and 
the confinement in z direction is modeled by a harmonic 
potential with frequency uj. In this gauge, the Hamilto- 
nian can be rewritten in the form 



H 



Px+Pz 

2m 



P'l 



2m 2 



(57) 



with the renormalized effective mass in y-direction, m = 
m(l + Aui 2 L I 'uj 2 ) , the renormalized confining energy Tiuj = 
hu>^/l + 4lu'^/uj 2 , and a shift in the confining potential 
zq = 2p y u! l i 'muj 1 which depends on the momentum p y in 
the y direction and the Larmor frequency uj, = eB/2mc. 
Note that the corrections due to the magnetic field in the 
resulting 2D Hamiltonian 



H 2D = 



Px 

2m 



Py_ 

2m 



(58) 



are of the order ui%/oj 2 or (a 2 /Zs) 4 , where Is = y/hc/eB 
denotes the magnetic length and a z — ym/mw the con- 
finement length. Usually, we are interested in the case of 
strong confinement and moderate magnetic fields where 
a 2 <C Ib, therefore m ss m up to small corrections. In 
this case, an in-plane magnetic field does not affect the 
orbital degrees of freedom of the 2D electrons. 

The condition (b) can be fulfilled in two ways. Either 
we have a theoretical description of the dependence of 
the Hamiltonian parameters (J, B,) on the control pa- 
rameters (V, a, B, E) or this relation is first mapped 
out in experiment and the obtained data is used later 
for the control of the device. A good approximate de- 
scription is possible in the case of adiabatic switching. 
In order to demonstrate this, we cast the microscopic 
Hamiltonian into the form H(t) = Hq + V(t). Then we 
find the instantaneous eigenstates |n(£)) and the corre- 
sponding instantaneous eigenvalues e n {t) by solving the 
time-independent Schrodinger equation for fixed time t. 
The instantaneous eigenstate |«(i)) is a good approxi- 
mation for the time evolution of the initial state |n(0)), 
provided the adiabaticity criterion 



(m\d t V\n) 



1 



(59) 



is met, where t s denotes the switching time. Eq. J5£ 
means that the change of the external control parame- 
ters during the switching time should be much smaller 
than the level spacing in the microscopic Hamiltonian. 
In the case of coupled quantum dots in the adiabatic 



regime, J(t) — e t (t) — e s (i) is the level spacing between 
the instantaneous singlet and triplet energies. 

Note also that if V(t) respects some symmetry, there 
can be selection rules that make Eq. (^9|) less stringent. 
In the case of two coupled quantum dots with an applied 
homogeneous magnetic field the total spin is conserved by 
V(t) and therefore only transitions to higher orbital lev- 
els of the quantum dots are relevant. Therefore, the less 
stringent condition l/r s ss | V'/V] <C Al/Ti is sufficient 



for adiabatic switching 16 1. Here, Ae denotes the orbital 
level distance averaged over the switching time. Since in 
this case the Zeeman energy is independent of the space 
coordinates, it commutes with the orbital Hamiltonian 
and does not affect adiabaticity. The case of inhomo- 
geneous magnetic fields is more intricate. The lack of 
a selection rule enforces the more stringent adiabaticity 
condition @ l/r s « \V/V\ < J/H < Ae/h, where J 
denotes the average exchange coupling during the switch- 
ing. In addition to this, the Zeeman term does not com- 
mute with the orbital Hamiltonian in the case of inho- 
mogeneous fields and therefore also influences J. Due to 
these difficulties, we presently do not know how to cal- 
culate the parameter J in Eq. (|l|) in the presence of an 
inhomogeneous field, Bi ^ B 2 . 

The condition of synchronous switching (c) is mainly 
a technical issue. We would like to stress that the choice 
of the pulse form has a decisive influence on whether the 
adiabaticity condition Eq. ( j59|) can be satisfied or not. 
It is quite easy to sec that a rectangular pulse is unsuit- 
able because it has infinite derivatives. Both Gaussian 
(exp(— t 2 / At 2 )) and exponential (exp(— |t|/A£)) pulses 
are far better than a rectangular pulse. The exponen- 
tial pulse has the advantage that |V"/V| is independent 
of t compared to the Gaussian pulse where | V"/V| cx t. 
However, the exponential pulse has the disadvantage that 
it has a cusp at t = which causes algebraically de- 
caying tails in its Fourier spectrum. We can combine 
the advantages of both pulses by using the sech pulse, 
sech(i/A£) = l/cosh(f/A£). Since all the pulses have to 
be cut off at some finite time ±t s /2, we choose the width 
of the pulse At smaller than the actual switching time r s , 
i.e. we choose a — t s / At > 1. By substituting the sech 
pulse into the adiabaticity condition Eq. (|59|), we obtain 
the condition r s 3> ah/Ae in the case where the spin is 
conserved (homogeneous magnetic field) and t s 3> ah/ J 
otherwise. 



VII. APPLICATIONS 

We will now give a detailed description of how a sys- 
tem of three coupled quantum dots could be controlled 
in order to test the functionality of three bit quantum 
error correction in that system. We denote the maximal 
coupling and magnetic field that can be applied by J max 
and -B max - If only one of the parameters Jy, B^ can 
be made non-zero at a given instant, then the following 
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serial-pulse sequence has to be applied: 



step 


duration 


parameter 


value 


1 




V 




2 


T//4 


Jl2 


J ni&x 


3 


tb/2 


Bl 

Z 




4 


T//4 

•J 1 


Jl2 




5 


TB /4 


Bl 

z 




6 


tb/4 


Bl 

z 


_^ maX 


7 


tb/4 


Bl 




8 


ts/4 


4 

V 




9 


T//4 


Jl3 




10 


r s /2 


Bl 


■^max 


11 


T.//4 


Jl3 


<Anax 


12 


rs/4 


Bl 


-^max 


13 


rs/4 


Bl 


^max 


14 


r s /4 


B 3 


^max 


15 


r n 


B x 


random 


16-29 


repeat 


1 - 14 





(60) 



where r,/ = 27r?i/J max and T£ = 2tt?i/ g[iBB max . Step 15 
describes the artificial introduction of noise into the sys- 
tem by applying a random magnetic field in the a; direc- 
tion, causing random spin flips in a time r„ < ir/gfisBx, 
where denotes the mean amplitude of the random B 
field. After step 29 is completed, qubits 2 and 3 are mea- 
sured and qubit 1 is flipped (by applying B\ = B max for 
time Tb/2) if both measurements yield 1 (spin down). 
The total switching time for steps 1 to 29 then amounts 
to T s = 6t b + 2tj + T n . 

In a device where parallel pulses are possible, i.e. where 
the conditions (a)-(c) from Sec. VI are fulfilled, the fol- 
lowing pulse sequence can be applied with the same ef- 
fect: 



i 


n 








Jl2 


B\ x 


B-2x 


Bly 


Biy 


Biz 


B 2z 


1 


.464 t 


.402 


-.054 


1 


.442 


.419 


-.905 


.851 


2 


M7t 


1 


-.356 


.559 


-.622 


.396 


.572 


.176 






Jl3 


B\ x 


B% x 


Bly 


B S y 


Biz 


Bzz 


3 


.464 t 


.402 


-.054 


1 


.442 


.419 


-.905 


.851 


4 


.617 t 


1 


-.356 


.559 


-.622 


.396 


.572 


.176 








B\ x 


B2x 


Bsx 








5 






rnd 


rnd 


rnd 








6 
1 


2.162r 


repeat 


1 

1 












9 






4 













(61) 

We have assumed that the maximal Zeeman energy is 
equal to the maximal coupling, g[J,BB niax — J max = M, 
and defined r = tb = Tj. All parameters are given 
in units of M. The parameters in every step can be 
multiplied by any pulse shape f(t) with /J" 1 f(t)dt = 1, 
where n denotes the duration of step i. Note that in 
every step, the pulse shape has to be the same for all 
parameters. Parameters that are omitted in Eq. ( |6l| ) are 
set to zero. The total switching time in this parallelized 



version amounts to r p = 4.3252 r + r n , compared to r s 
8t + r n in the case of serial switching. 



VIII. CONCLUSION 

We have studied the minimal requirements for the im- 
plementation of the XOR gate, the conditional phase 
flip (CPF) gate, the encoding circuit E used for three 
bit error correction, and the teleportation encoder Et, 
all for Heisenberg-coupled spins with s = 1/2. In ad- 
dition to this, we have also considered anisotropic spin- 
spin coupling as described in the XY model. Two differ- 
ent methods for generating quantum gates with a time- 
dependent Hamiltonian have been discussed and com- 
pared, the "conventional" serial pulse method and a new 
method involving parallel pulses. 

The main results of our work are the parallel pulses for 
the conditional phase flip (Eq. (§f)) and XOR (Eq. @) 
using Heisenberg dynamics, and the corresponding re- 
sults (Eq. (||) and Eq. (|J)) for XY dynamics. The 
direct parallel-pulse sequence Eq. (^) for the three-bit 
encoder E was found; however, it is possible that a faster 
pulse sequence for this gate can be found with more nu- 
merical effort. 

The following results for serial switching have been 
found: There is an analogue of the known circuit Eq. (j^) 
for CPF (cf. Fig. 0) for systems with XY coupling, which 
is given in Eq. (Eq). For Heisenberg coupling, we have 
proved that the known circuit Eq. (o) is optimal in the 
sense that CPF cannot be obtained with one 'square- 
root-of-swap' gate. For the proof we invoked the set V(U) 
of all product states that are mapped onto product states 
by a quantum gate U ; V(U) helps to distinguish quantum 
gates modulo concatenation of single-qubit gates. The 
same tool was also used to prove that the encoder E for 
quantum error correction cannot be generated with se- 
rial switching with fewer than four 'square-root-of-swap' 
gates. The same is true for the encoder Et for the tele- 
portation of one qubit. 

The results for the parallel-pulse XOR for isotropic 
Heisenberg interactions and the results for CPF and 
XOR for XY interactions, Eqs. @s ©, and (g), 
and for the three bit encoder Eq. (Q) were all found 



using the computer algorithm described in Section |IV. 
This algorithm searches for a (parallel) pulse sequence 
for an arbitrary quantum gate operating on any number 
of qubits. The number of qubits and the complexity of 
the pulse sequence that can be studied are only limited 
by the available computational resources. 

Quantum computations are very often presented in the 
form of quantum circuits, i.e. as a sequence of gates be- 
longing to a small set of universal gates. Our examples 
of parallel-pulse gates illustrate that such quantum cir- 
cuits are in general not the most efficient way of per- 
forming a quantum computation. The reason for this is 
that standard quantum circuits only allow the use of a 
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small fraction of the possible time evolutions that can 
be generated by the underlying Hamiltonian. While for 
the two- and three-qubit gates we have studied here, we 
could optimize the switching time by typically a factor of 
about two by using parallel pulses, it can be speculated 
that for gates operating on many qubits or whole quan- 
tum computations, switching times could be reduced by 
a much larger amount. Note also that the parallel pulses 
we have studied here represent only a small subset of the 
possible time evolutions themselves, since we have been 
restricted to very simple discrctized pulses of up to three 
time-steps. 

While quantum circuits are very intuitive and provide 
an excellent framework for the theoretical study of quan- 
tum algorithms and their connection to classical algo- 
rithms, the representation of quantum gates or whole 
computations as parallel pulse sequences may turn out 
to be more efficient for a number of physical implemen- 
tations. 
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